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^ ' Abstract. The pre-main sequence binary system V773 Tau A shows remarkable flaring activity around periastron passage. 

^^ I Here, we present the observation of such a flare at a wavelength of 3 mm (90 GHz) performed with the Plateau de Bure 

Cn ■ Interferometer.* We examine different possible causes for the energy losses responsible for the e-folding time of 2.31± 0.19 

^~|^ ' hours of that flare. We exclude synchrotron, coUisional, and inverse Compton losses because they are not consistent with 

JN^ , observational constraints, and we propose that the fading of the emission is due to the leakage of electrons themselves at each 

reflection between the two mirror points of the magnetic structure partially trapping them. The magnetic structure compatible 

with both our leakage model and previous observations is that of a helmet streamer that, as in the solar case, can occur at the top 

of the X-ray-emitting, stellar-sized coronal loops of one of the stars. The streamer may extend up to ~ 20 R, and interact with 

the corona of the other star at periastron passage, causing recurring flares. The inferred magnetic field strength at the two mirror 

M-(| points of the helmet streamer is in the range 0.12 - 125 G, and the corresponding Lorentz factor, y, of the partially trapped 

r-^ . electrons is in the range 20 < y < 632. We therefore rule out that the emission could be of gyro-synchrotron nature: the derived 

$^ ' high Lorentz factor proves that the nature of the emission at 90 GHz from this pre-main binary system is synchrotron radiation. 

C^ ■ Key words, stars: coronae, stars: individual: V773 Tau, stars: flare, stars: pre-main sequence. Radio continuum: stars 
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1. Introduction the magnetic field intensity and the electron Lorentz factor can 

^ be derived, yielding information on the emission mechanism. 

. . . The dynamo theory (Parker 1955) explains how differential ro- Concerning the first issue, a, one can distinguish among 

tation generates a toroidal field in the interior of a star from ^j^^-gg possible origins for flares- 
an initial stellar dipole field, and how convection, bringing this 

field up to the surface, creates the coronal magnetic arc-like 1- In a i/ng/e star, flares occur when new emerging loops from 
structures called loops. A close relationship between magnetic below the stellar surface move into older, already-existing 
loops and flares exists. In fact, the physical mechanism invoked loops (Heyvaerts et al. 1977; Massi et al. 2005). Such in- 
to explain solar and stellar flares today is magnetic reconnec- truding loops have indeed been observed on the Sun in 
tion, which occurs when different loops interact with each other high-resolution maps obtained with the Nobeyama Radio 
or when field lines of the same loop are stretched. The re- Interferometer (Nishio et al. 1996). 

leased magnetic energy accelerates a part of the thermal elec- 2. In protostars or young (pre-main sequence) stellar systems, 

irons trapped in the loops to high energies, causing the flare which are stiU surrounded by a dense accretion disk, flares 

(Golub & Pasachoff" 1997; Priest & Forbes 2002). can occur because of stretching, disruption, and reconnec- 

We have observed a strong flare at mm wavelengths toward tion of magnetic field lines between the star and its disk 

the T Tauri binary system V773 Tau A. In the context of the (Feigelson & Montmerle 1999). 

outlined scenario, our aim is to investigate: fl) the flare location 3. In close binary systems consisting of late spectral type 

and hence the geometrical structure of the involved magnetic stars, much higher magnetic activity is expected compared 

field, and b) the nature of the observed emission, i.e. whether to the Sun (because tidal synchronism increases the rota- 
tion rate and therefore the dynamo efficiency). Moreover, 
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(Germany), and IGN (Spain). the young binary system V773 Tau A, which shows a high 
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rate of radio flares around the periastron passage (Massi et 
al. 2002). 

Concerning the second issue (b), the usuaUy observed radio 
emission from steUar coronae shows spectra peaking around 
10 GHz and is consistent with the gyro-synchrotron process 
(Mutel et al. 1985; Bastian et al. 1998; Giidel 2002). However, 
some solar and stellar flares can also be observed at millimeter 
wavelengths (Fig. 1), and the spectra are thought to be due to 
the superposition of a gyro- synchrotron spectrum with a syn- 
chrotron spectrum, the last one peaking in the shorter sub- 
millimeter to the far-infrared range (Kaufmann et al. 1988; 
Kaufmann et al. 2004). Gyro-synchrotron emission, associated 
with mildly relativistic particles, i.e., with a Lorentz factor of 
■y < 10, is circularly polarized. In contrast, synchrotron emis- 
sion comes from highly relativistic particles (y » 1) and is 
linearly polarized (Dulk 1985; Phillips et al. 1996; Tsuboi et 
al. 1998). 

To probe the existence of relativistic electrons in the V773 
Tau A system and to constrain their energy and the strength 
and topology of the magnetic field, we observed it at millimeter 
wavelengths using the IRAM Plateau de Bure Interferometer 

The paper is organized as follows. I^rcviews previous re- 
sults on V773 Tau A whereas in |3lour new observations and 
results are presented. In l^we evaluate which kind of energetic 
processes may cause the trapped electrons to lose their energy 
and derive the Lorentz factor and magnetic field strength. These 
deductions are confirmed in ^by the analysis of the flare ris- 
ing time, ^discusses the magnetic field configuration, while 
^presents our conclusions. 
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Fig. 1. Examples of spectra of the Sun (squares: Akabane et al. 
1973; empty triangles: Zirin and Tanaka 1973) and UX Arietis 
(filled triangles: Beasley & Bastian 1998) exhibiting a "flatten- 
ing" towards mm- wavelengths. Notice that the vertical scale is 
in both solar flux units (for the solar spectra) and mJy (for UX 
Arietis). 



2. The V773 Tau A binary system 

V773 Tau, at a distance of 148±5 pc (Lestrade et al. 1999), is 
a quadruplet of T Tauri stars within an area of radius less than 
100 AU (Woitas 2003). The system has two stars with orbital 
periods of years and two shorter period inner stars. The object 
of our studies is the V773 Tau A binary, with an orbital period 
of only 51.075 days (Welty 1995). The orbit of V773 Tau A is 
moderately eccentric (e=0.3) and the epoch of periastron pas- 
sage, corresponding to orbital phase O = 0, is r()=2449330.94 
JD. The stars, both very active, have similar radii (/?,= 2.4 Rq) 
and rotation periods (about 3 days), but seem to have different 
masses (see Welty 1995). 

The system possesses only weak hydrogen emission lines 
in its spectrum, which implies that its components are weak- 
line T Tauri stars with mass loss. This follows from the pres- 
ence of forbidden emission lines (Cabrit et al. 1990). The pres- 
ence of a disk is controversial (Skinner et al. 1997) since the 
observed infrared excess may also be attributable to V773 Tau 
D, the forth member of the quadruplet (Duchene et al. 2003). 
Also, it is difficult to establish any thermal component in the 
millimeter emission (which would imply a disk), since it is 
strongly dominated by a variable non-thermal component typi- 
cal of flaring coronal activity (Skinner et al. 1997). At 2.7 mm, 
variability was found with the flux density changing from 30 
mJy to less than 3 mJy in a period of a few months (see Skinner 
et al. 1997). Beckwith et al. (1990) measured a flux density at 
1.3 mm, S =42 + 6 mJy, whereas an observation several years 
later gave S =24 + 4 mJy (see Skinner et al. 1997). With 
the Eff'elsberg 100-m telescope, Massi, Menten, & Neidhofer 
(2002) observed a flare at 7 mm with S = 68 ± 13 mJy, and 
three months later they only obtained an upper limit (3cr) of 23 
mJy. 

V773 Tau A is not the only pre-main sequence object with 
mm-flaring emission. For example. Bower et al. (2003) ob- 
served a giant outburst at 86 GHz in the Orion source GMR-A, 
also a weak-line T Tauri star. Long -lasting (13 days) millime- 
ter activity in that source was observed by Furuya et al. (2003). 
However, V773 Tau A is the only known system with coronal 
solar-like magnetic activity where linearly polarized emission 
together with the more common circularly polarized emission 
was measured (Phillips et al. 1996). The high level of variable 
millimeter emission together with the presence of linear po- 
larization indicates that a quite energetic electron population is 
present. The fact that these energetic electrons have rather short 
lifetimes is worthy of note for the following analysis. Phillips et 
al. (1996) find that total intensity variations at different wave- 
lengths and epochs share common time scales and establish a 
recurring interval for the decay of radio events with a typical 
1/e time scale of only about 1 hour 

The two components of V773 Tau A are both active (Welty 
1995); flares on both of them are expected due to mechanism 
a-l discussed in the introduction. In addition, however, there 
is evidence for mechanism a-3 (inter-binary loop collision): by 
folding the data, monitored over 522 days, with the orbital pe- 
riod (5 1 .075 days), the flares appear to be clustered at perias- 
tron passage (Massi et al. 2002). Around the periastron, more 
than one interbinary interaction seems to occur: three consec- 
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utive flares, separated by a time interval of 3 - 4 days, have 
been observed during a continuous observation around perias- 
tron passage (Fig. 3 in Massi et al. 2002). The hypothesis of a 
stable active region simply eclipsed by the body of the star dur- 
ing each rotation (Pmt= 3 days) is ruled out by the short life- 
times (1-2 hour only) of the energetic electrons (as discussed 
above and in §3). On the contrary, a large asymmetric structure 
rooted on one of the two rotating stars could explain the oc- 
currence of the three consecutive flares, as repeated collisions 
(at each rotation) with the corona of the other star. The extent, 
H, of this magnetic structure must be large enough to allow 
the observed consecutive collisions around periastron passage. 
The periastron separation is 56 solar radii, Rq, which for a stel- 
lar radius, /?. = 2.4 R^ (Welty 1995), corresponds to 23 R,. 
Indeed, a large magnetic structure has been imaged by using 
Very Long Baseline Interferometry (VLBI, Phillips et al. 1996; 
orbital phase - 0.1) at /I - 3.6 cm, with two peaks of dif- 
ferent intensity separated by 0.17 AU (15 /?,). Observational 
constraints for the size of the emitting region are therefore in 
the range 15-20/?,. 

The fact that linear polarization has been observed implies 
that Faraday depolarization is not effective and puts an upper 
limit to the density (n < 10^ cm"^) of the plasma confined in 
this extended magnetic structure (Phillips et al. 1996). X-ray 
observations give evidence not only of a density much above 
this limit, but also of the existence of smaller, solar-like, coro- 
nal loops in V773 Tau A. Skinner et al. (1997) interpreted the 
light curve of a hard X-ray flare in V773 Tau A as being due to 
the rotational modulation of the emitting flaring region, deter- 
mining a size of H < 0.6 R». Tsuboi et al. (1998) interpreted the 
decay of another hard X-ray flare as being due to radiative cool- 
ing, obtaining a size of 1 .4 /?«. Therefore, radio and X-ray emis- 
sion come from spatially separated regions: denser and smaller 
ones associated with X-ray emission and larger and more dif- 
fuse ones associated with radio emission. Further confirmation 
of the presence of two structures comes from the multiwave- 
length campaign on V773 Tau A carried out by Feigelson et 
al. (1994) showing radio variability combined with a steady X- 
ray flux. Using our millimeter observations, the following sec- 
tions are aimed to better define what kind of extended magnetic 
structure might confine the radio-emitting plasma, the intensity 
of the magnetic field, and the energy of the trapped particles. 

3. Observations and data reduction 

V773 Tau A was observed with the IRAM Plateau de Bure 
Interferometer at 1 mm and 3 mm around the two periastron 
passages in August and November 2003. The configuration 
was 5D (five antennas) in August and 6Cp (the six-antenna 
configuration) in November Observations were done in LCP 
for the 3 mm observations and in RCP for the 1 mm observa- 
tions. The continuum was measured with four correlator units 
per receiver, twice with a bandwidth of 320 MHz and twice 
with a bandwidth of 160 MHz. Cyclical observations (every 
20 minutes) of the sources 0528-1-134 and 3C286 were used 
for phase and amplitude calibrations. The data were reduced 
with the programs CLIC and MAPPING of the GILDAS soft- 
ware package developed by the Grenoble Astrophysics Group. 



Table 1. Flux density in November 2003 



Day 



, (mjy) 



, (mJy) 



13 
14 



5.0: 
6.0: 



1.2 
1.2 



3.0: 
5.0: 



0.4 
0.3 



The source was always unresolved. A complete flare (Fig. 2) 
was observed on August 6 (O = 0.1). With about 12 mm of 
precipitable water vapor in the atmosphere, no observations at 
A - I mm were possible, and those at A - 3 mm had the high 
rms noise level of 17 mJy; the 3 mm receiver was tuned to 90 
GHz and had a system temperature around 1000 K. The best 
power-law fit to the decay part of the light curve shown in Fig. 
2 gives an e-folding time of t = 2.31 ± 0.19 hour Around the 
periastron passage in November, we were able to observe un- 
der very good weather conditions for two consecutive days at 
<[) = 0. The flux density was at a quiescent constant level with 
a nearly flat spectrum (Table 1). 



300 - 





Time [hours] 



Fig. 2. V773 Tau A flare observed August 6, 2003 (O = 0.1) 
with the Plateau de Bure Interferometer at /I = 3 mm. The error 
bar corresponds to 2cr. 



4. Rapid flare decay: leakage of the emitting 
particles 

As outlined in the introduction, the emission observed at 
3.3 mm (90 GHz) can be interpreted as synchrotron radia- 
tion. In synchrotron radiation, the emission of each relativis- 
tic electron with Lorentz factor y moving in a magnetic field 
B (in Gauss) is centered around a peak spectral frequency, vq 
(Ginzburg and Syrovatski 1965) of 



vo = 1.8xlO^By2 



Hz. 



(1) 



This implies that, to reproduce the observed emission at v = 90 
GHz, electrons must exist for which the following relationship 
holds: 



By^ = 5 X 10'*. 



(2) 
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The temporal evolution of the flux density shown in Fig. 2 
clearly indicates that the emitting electron distribution, respon- 
sible for the observed synchrotron emission, is subject to some 
losses. In particular, the derived e-folding time of 2.3 1 hours for 
the mm flare is of the same order of that measured by Phillips 
et al. (1996) for cm emission (about one hour). This time scale 
implies that the electrons must be at least partially trapped in 
the emitting region, which would otherwise be rapidly depleted 
in a time equal to the size of the region divided by the speed of 
the electrons: 20/?,/c ~ 112 s. It is well-established that par- 
ticles trapped in magnetic structures spiral around field lines 
to a point (the "miiTor point") where the magnetic field lines 
converge and the field strength is sufficient to cause the par- 
ticles to reverse direction and travel back to the other mirror 
point. Reflected back and forth between mirror points (known 
as "bounce" motion), the particles remain trapped and, contin- 
uously spiraling, emit synchrotron radiation until they have lost 
their energy (Roederer 1970). 

However, the decay time of 2.31 hours is not attributable 
to synchrotron losses. As derived in the Appendix, synchrotron 
losses would explain an intensity decay time of that duration 
only for an emitting region of size < 3.8 R«, whereas the ob- 
served size (§2) is H > l5Rt. Neither inverse Compton losses, 
due to the radiation field of the star, and/or to the synchrotron 
emission of the flare (synchrotron self-Compton (SSC)), can be 
responsible for the observed decay; Compton losses are even 
lower than synchrotron losses (see Appendix). On the other 
hand, even collisional losses cannot account for the observed 
decay time, since collisions are important only for high density 
plasmas (see Appendix) in contrast to the upper limit on the 
plasma density « < 10^ cm"^ derived from Faraday depolariza- 
tion (§2). 

Since the fast decay of the emission cannot be attributed to 
energetic losses of the electrons, it may be caused by leakage 
of the particles themselves, meaning that at each bounce mo- 
tion at one mirror point, some of the electrons are able to leave 
the trapping region and escape into free space. This is possible 
for electrons having a pitch angle, (the angle between their 
velocity direction and the magnetic field vector), smaller than 
the loss cone angle 



sm 6/ < sm ( 



HI) 



1/2 



(3) 



Here, Bi is the magnetic field intensity far out in the corona 
at a distance (Hi), while Bq is that at a lower distance (Ho). 
Pitch angle scattering, also known as "pitch angle diffusion" 
(Mekose & Brown 1976), is one of the main processes by 
which trapped particles are lost. Weak and strong diffusion 
have been introduced by Melrose & Brown (1976). However, 
in order to be effective, weak diffusion demands a high-density 
medium (« > lO'^cm"-'), in contrast to what we discussed 
above. Also, Lee et al. (2002) found that for the microwave 
decay of a solar flare, the assumption of strong diffusion seems 
to be the most appropriate. In this scenario, the emission de- 
creases by a factor of 1/e if the precipitation rate 



1 



2 °(i/i-//n) 



(4) 



(Melrose & Brown 1976; Lee et al. 2002) times the decay time 
(2.31 hours, in our case) equals unity. This condition, valid for 
small Oq values, implies: 



742 X 



Bi 



R* 



Bo Hi - H 



1. 







(5) 



For a magnetic dipole field, Bi/Bo = [Hq/Hi]^, possible solu- 
tions of Eq. Q are values in the ranges Ho ~ (2 - 5)Rt and 
Hi ~ (10 - 20)Rt. Assuming a magnetic field strength at the 
stellar surface of the order of 1 kG (see Appendix), the corre- 
sponding values of the magnetic field intensity are in the range 
125 - 8 G for Hq and 1 - 0. 12 G for Hi . This scenario is consis- 
tent with synchrotron emission, since from Eq. (|3, the corre- 
sponding Lorentz factor of relativistic electrons emitting at 90 
GHz is 20< r < 632. 

This scenario is also in accordance with the VLBI ob- 
servation of Phillips et al. (1996), since the analogon of Eq. 
Q derived for the decay time of one hour and solved for 
Hi - Ho - 15/?. with Hi ~ 20/?, gives B in the range of 
4.6 - 0.12 G. These values of the magnetic field for v - 8.2 
GHz (Phiflips et al. 1996) coiTespond (see Eq. (1)) to 31 < y < 
195 indicating, therefore, the presence of relativistic electrons 
in agreement with the observed linear polarization. 

5. The flaring phase: a propagating shock 

In the previous section, we established the most likely physical 
process responsible for the rapid flare decay. The aim of this 
section is to test whether the additional constraint set by the 
flare rising time confirms the previously obtained results. 

Recent solar flare studies (Tanuma & Shibata 2005; Asai 
et al. 2004) have shown that the observed downflow motions 
can create fast shocks that are related to observed non-thermal 
bursts, since shocks are supposed to be sites of electron acceler- 
ation (Aschwanden 2002). In our scenario, in which reconnec- 
tion takes place far out where the two stars' magnetospheres 
interact (i.e. at Hi), the shock induced by magnetic reconnec- 
tion can propagate down to Hq along the magnetic structure at 
the local Alfven velocity 



va(H) = 



B(H) 



, -7 10''l — 

y[AmnHn{H) \R 



H 



cm/ sec 



(6) 



(for a dipole magnetic field with B(/?,) = 1000 G and a particle 
density n ~ 10^(////?,)"^ cm"-'), thus inducing successive par- 
ticle acceleration events. The shock propagation speed can be 
related to the flare rise time, fiise. 



mse 



•JHo 



dH 



(7) 



from which the following relationship can be derived: 



Hi? ^ /^^ lBiR^\ /^^\- //^\^ 
rJ \hour/\103G/\10''cm-V \rJ 

For frise, - 4.5 hours and Ho between (2 - 5) /?,, the resulting 
height Hi is in the interval of (13 - 27) /?,, assuming that B(R,) 
is in the range 1000-3000 G (see Appendix) and the density 
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Fig. 3. Top: Variations of the height Hi (Eq. (8)), where mag- 
netic reconnection occurs, as a function of magnetic field inten- 
sity and plasma density (see ^. Bottom: Alfven radius (Eq. 
(9)) as a function of magnetic field intensity and mass loss rate 
(see© 



is in the range 10 - 10 cm~ , as in Fig. 3. This range for 
Hi is well-consistent with the scenario of a reconnection event 
at a large stellar distance and in agreement with the resulting 
interval (10 - 20)Rt from the leakage model. 



6. Magnetic configuration 

In the previous sections, in the context of our leakage model, 
we established magnetic field values in the range 125 - 8 G (for 
Ho ^2-5 R,) and 0.12 G (at //j = 20 R,) at the two mirror 
points of the confining magnetic structure. Also, we determined 
values of 4.6 G (at Hq ^ 6 R,) and 0. 12 G (at Hi = 20 R,) for 
the centimeter observations of Philipps et al. (1996). In both 
cases, in a high-resolution image of the source brightness dis- 
tribution, the two miiTor points would coiTespond to two peaks 
(displaced by Hi - Hq) of rather diff'erent intensity, because of 
the relationship between emissivity and B (Dulk 1985). This 
predicted morphology finds its confirmation in the Very Long 
Baseline Array (VLBA) image of Phillips et al. (1996). The 
VLBA image, given in the small box of Fig. 3, shows a strong 
brightness distribution peak, which we may identify with the 



first mirror point, Pq, separated by 15 R* from a second weaker 
peak, which may be coincident with the second mirror point. 

Pi. 

Let us assume the scenario of two stellar coronae inter- 
acting at periastron passage. One of the two coronae must be 
very asymmetric, in accordance with the asymmetric magnetic 
field configuration implied by the Massi et al. (2002) obser- 
vations of three consecutive flares around periastron (see Sect. 
2). On the other hand, a single giant stellar loop is ruled out by 
X-ray observations which, on the contrary (see Sect. 2), pro- 
vide evidence for a stellar-sized coronal loop (Skinner et al. 
1997; Tsuboi et al. 1998). However, in the case of the Sun, ex- 
tended magnetic structures are sometimes present above coro- 
nal loops. Indeed, when the Sun is observed with white-light 
coronographs, one can see that above the top of some coro- 
nal loops (Fig. 4), a dome-shaped structure (the "helmet") is 
formed, extending out to 2 - 4 solar radii, with streamers that 
can extend out to many solar radii (Suess & Nerney 2004; 
Endeve et al. 2004). A helmet streamer is therefore a plausible 
physical model for the elongated asymmetric structure implied 
by the previous observations of Massi et al. (2002). Moreover, 
a helmet streamer structure with its two miiTor points at Pq and 
Pi (Fig. 4) matches the results of our leakage model based on 
PdBI data, as well as the appearence of the source in the VLBA 
image. 

Finally, we point out that the helmet streamer in the sim- 
plified sketch of Fig. 4 is represented, for simplicity's sake, by 
a straight structure. In real highly conducting stellar plasmas, 
where the ionized gas is frozen to the field lines, more complex 
morphologies will occur: wherever magnetic pressure exceeds 
gas pressure, the fluid moves only along the field lines. In con- 
trast, when gas pressure dominates, the magnetic field lines do 
move or bend, depending on the fluid motion. In a rotating star 
with mass loss rate M and terminal wind velocity v^, the dis- 
tance (called Alfven radius, Ra) where the corotation of the 
wind can no longer be enforced by the magnetic field and the 
magnetic field lines become curved is, as given by Andre et al. 
(1988), 



R 



^=26(^V 
R, UO^G/ 



M 



10-iOMo yr- 



{S)'« 



.lO^kms"', 

Adopting a range lO^'^-lO^'^M© yr"' for the mass loss rate of 
a weak T Tauri star (Andre et al. 1992), and a wind terminal 
velocity of 300 km/sec (Ardila et al. 2002), we can conclude 
that the Alfven radius for P = 3 d (Fig. 3 bottom) is in the 
rangeof(14-30)/?,. 

7. Conclusions 

Our previous observations of the binary system V773 Tau A 
(Massi et al. 2002) have shown inter-binary loop collisions at 
periastron passage. We performed millimeter wavelength ob- 
servations aimed to better define the extended magnetic struc- 
ture confining the radio-emitting plasma, its topology, strength, 
and the energy of particles trapped in it. Based on our results, 
we propose the following scenario. 

Radio and X-ray emission from V773 Tau A arise from spa- 
tially separated volumes: X-ray emission arises from denser 
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Fig. 4. Sketch of the magnetic structure of a helmet streamer, 
above a coronal loop, as treated in magnetohydrodynamic mod- 
els (Endeve et al. 2004). The small box contains the VLBA im- 
age by Phillips et al. (1996). 




at apoastron (40 /?,). The size of the helmet streamer should 
therefore be comparable with the periastron distance, but still 
be less than 40 /?,. In this paper, where we report a flare at 3 
mm with 5peak = 360 ± 17 mJy, we have derived the size "Hi" 
of the emitting region (i.e. the helmet streamer) in three com- 
pletely independent ways. Interpreting the rising time of the 
flare in terms of propagation towards the low corona of a shock 
due to a magnetic reconnection event occurring at large height, 
we have determined a height Hi in the range of (13 - 27) /?,. 
This size agrees with that determined for the Alfven radius, i.e. 
in the range (14 - 30) /?,. By interpreting the rapid decay of 
the flare as leakage of the emitting electrons from the trapping 
magnetic structure of size Hi, we determine Hi in the range 
(10 - 20) /?«. As a result, one of the two stars of the system 
might have a corona of a few stellar radii only, and interbinary 
collisions can still occur because of a large (i.e. average size 
Hi ~ 19^,) helmet streamer located on the other star. 

The leakage model has shed light on the nature of the elec- 
trons responsible for the flare. We state that interbinary colli- 
sions are able to accelerate electrons to relativistic velocities. 
In fact, we determine a Lorentz factor y up to 632. This con- 
firms that the emission in this pre-main sequence system is syn- 
chrotron radiation, as postulated from the observed linear po- 
larization (Phiflips et al. 1996). 

In particular, streaming back and forth between the two 
mirror points (Pq and Pi in Fig. 4) of the helmet streamer, 
with magnetic fields between 0.12 G and 125 G, the relativis- 
tic electrons produce the observed synchrotron radiation. The 
trapping by the helmet streamer is, however, only partial and 
during each reflection particles escape and the emission fades 
out relatively quickly. Leakage explains the puzzling finding 
that flares have similar decay times (1-2 hours), not only in 
different radio bands, but also in the millimeter band. The latter 
similarity would be difficult to explain on the basis of any loss 
mechanism dependent on the electron energy. 



Fig. 5. Sketch of the binary system with the orbit parameters 
of Welty (1995). One star with helmet streamers is shown at 
different positions along the orbit. Interbinary collisions may 
occur at periastron passage (see text). 



and smaller regions (size ^ 1/?,) like the closed loop region 
under Pq in the sketch of Fig. 4 (see §2). The radio emission 
emerges from larger and more diffuse regions that are similar 
to the solar helmet streamers. A helmet streamer from one of 
the stars, formed at the top of a coronal loop, as in the case of 
the Sun, could interact with the corona of the other star and 
could produce the strong radio flares observed around peri- 
astron passage (Massi et al. 2002). In this scenario, the lack 
of flares at some periastron passage (as in the PdBI observa- 
tions of Table 1) could be explained as the absence of a col- 
lision, possibly because no large enough helmet streamer was 
emerging on the hemisphere facing the other star. As shown in 
Fig. 5, because of the slightly eccentric orbit (e ^ 0.3, Welty 
1995) the distance at periastron (23 /?,) is appreciably less than 
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Appendix A: energetic considerations 

Accelerated non-thermal electrons are subject to energy losses 
both because synchrotron radiation is emitted and because of 
collisions with thermal electrons. In addition, energy losses are 
also due to the interaction of non-thermal electrons by inverse 
Compton (IC) processes with the stellar radiation field (exter- 
nal IC) and even with photons emitted by the synchrotron pro- 
cess itself in the radio band (synchrotron self-Compton (SSC)). 
If synchrotron losses are responsible for the observed de- 
cay (t, in hours), the following relationship holds between 
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magnetic field intensity B (in Gauss) and Lorentz factor y 
(Blumenthal & Gould 1970): 

7 2.2 xlO^ 
yB^ = . (A.l) 

^synchrotron 

From Eqs. jA.ll l and Q, with Tsynchrotron - 2.31 hours, the 
value of the magnetic field strength in the emitting region and 
the value of the Lorentz factor of the emitting electrons can be 
derived as: B = 56 G and 7=30. 

A magnetic dipole field, with B(Rt) at the stellar surface, at- 
tains the intensity of 56 G at a distance H/Rt - [B(Rt) / 56]' . In 
the literature, there are very few estimates of the magnetic flux 
B{Rt) in weak T Tauri stars. Basri et al. (1992) determined the 
upper limit of 1500 G for the product B/, where/ <1, for TAP 
35. An even higher value (B/ > 2000 G) is derived by Guenther 
et al. (1999) in the source LkCa 16. Basri et al. (1992) exclude 
fields above 3000 G, and more recently (Bower et al. 2003), the 
value of B=2600 G has been measured in a weak T Tauri in the 
Orion Nebula. Therefore, 1000 - 3000 G seems to be a likely 
range. For B{Rt) < 3000 G, the magnetic dipole field attains the 
intensity of 56 G at a distance H/Rt < 3.8 from the stellar sur- 
face. This scenario cannot account for the observed extended 
radio emission (§2). We also exclude synchrotron losses as be- 
ing responsible for the lifetime t - I hour of the Phillips et al. 
(1996) observations at vq = 8.2 GHz. The appUcation of Eqs. 
Q and ( lA.U would give B ~ 223 G and j - \, and hence 
emission coming from a region relatively close to the stellar 
surface (H/R^ < 2.4) and from mildly relativistic electrons, in 
contrast with previously reported considerations (§1 and 2). 

IC losses cannot account for the observed emission decay 
either. In fact, as derived hereafter, IC losses are even lower 
than synchrotron losses because the radiation energy density 
(U) is smaller than the magnetic energy density i.e., 



U < B^/9,7T. 



(A.l) 



A dipole magnetic field, with the intensity at the surface of the 
star in the range quoted above (say 2 lO-' G), has B^/^n ^ 1 .6 x 
lOHH/R^y^ erg cm-^ while U = LJ4kcH^ = 0.86 (H/R,)-^ 
erg cm"^ for known values of radius and luminosity of V773 
Tau A (R, ~ 2ARq, U ~ 3Lq) (Skinner et al. 1997; Welty 
1995). Hence, if the radiation of the star is the source of the 
seed photons (IC), condition JA.2> holds for H < 21/?». On 
the other hand, for SSC, the related flare luminosity is Lflaie = 
6.4 X 10^^ erg sec"', which, supposing that the flare takes place 
in a spherical region of radius aR^, with a > 0.01, gives U = 
6 10"^a"^ erg cm"^. Hence, condition ( IA.2> holds throughout 
the region of interest, since (H/R^) < 118 a'^^ is verified for 
every (H/R,) < 25. 

Finally, the thermalization time due to collisions of rela- 
tivistic electrons in an ionized gas of density n is given by 
(Petrosian 1985; Massi & Chiuderi-Drago 1992): 



Tcollisions - 4.16 X 10 
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(A.3) 



where n is in cm"^ and t is expressed in hours. Therefore, colli- 
sions would be important only for n~7'1.8xlO^ cm"^, which, 
for the relativistic electrons (y >>1) implied by the linear po- 
larization, does not match the condition derived from Faraday 
depolarization (n < 1.x 10^ cm"-'; see §2). 
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